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We investigate clogging in depth filtration, in which a dirty fluid is "cleaned" by the trapping of 
dirt particles within the pore space during flow through a porous medium. This leads to a gradient 
percolation process which exhibits a power law distribution for the density of trapped particles 
at downstream distance x from the input. To achieve a non-pathological clogging (percolation) 
threshold, the system length L should scale no faster than a power of lnui, where w is the width. 
Non-trivial behavior for the permeability arises only in this extreme anisotropic geometry. 
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Depth filtration is a mechanism for separating sus- 
pended particles from a carrier fluid by passing it through 
a porous medium [jl] [|] . The medium promotes efficient 
filtering both by increasing the area available for deposi- 
tion of suspended particles, as well as the exposure time 
of the suspension to the absorbing interfaces. This mech- 
anism is therefore widely used in a variety of biologi- 
cal, chemical, and engineering separation processes [yj. 
Depth filtration also raises basic issues in porous me- 
dia transport, as the medium becomes progressively con- 
stricted by particle trapping events so that a steady state 
is not achieved. This feedback between flow and structure 
governs the essential properties of filtration. Each pore 
blockage results in a small reduction of the permeability 
of the medium and ultimately a clogging (percolation) 
threshold is reached where the filter permeability van- 
ishes. Previous studies indicated that the permeability 
vanishes as a power law near the threshold, with an ex- 
ponent different from that of classical percolation 
From a practical perspective, filter performance is de- 
graded by particle trapping, so that understanding this 
trapping rate is of paramount importance. 

In this Letter, we formulate a geometrical description 
for depth filtration which provides intuition for the clog- 
ging process and leads to phenomenology outside the 
realm of classical percolation. For overlapping distri- 
butions of pore and particle radii, the trapped particle 
density distribution varies as a power law in longitudinal 
coordinate. Such a distribution should be readily ob- 
servable, for example, when opaque particles are passed 
through a glass bead-pack porous medium || . This gra- 
dient leads to an unusual percolation process where the 
percolation threshold (fraction of open pores) is very 
close to unity and where the permeability does not have 
power law behavior. 

Of the many microscopic interactions that underlie fil- 
tration, we focus on size exclusion where a par- 
ticle of radius r par ti c ie is trapped within the first pore 
encountered whose radius satisfies r porc < r par ti c ie- This 
size exclusion is the dominant effect in processes such as 
gel permeation in porous media and liquid chromatogra- 
phy. While other influences, such as Van der Waals forces 
between particles and pore surfaces, hydrodynamic and 



electrostatic interactions, etc., may be important, their 
faithful modeling is complex |[| and makes it difficult to 
identify the governing mechanism for a given macroscopic 
property. Our approach is to retain size exclusion as the 
only trapping mechanism in a geometric modeling of fil- 
tration and develop physical intuition for clogging from 
this idealized description. 
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FIG. 1. The bubble model. Each line represents fluid car- 
rying pore. 

We first introduce the quasi one-dimensional "bubble" 
model to describe filtration. This system consists of L 
links in series, in which each link is a parallel bundle 
of w bonds and each bond represents a pore (Fig. 1). 
This model can be viewed as a square lattice in which 
all perpendicular bonds are "shorted" . This system was 
introduced to account for the breaking of fibers (^] and 
to determine extremal voltages associated with break- 
down in resistor-fuse networks ^| . An appealing feature 
of this model in the context of percolation is that it ex- 
hibits finite-dimensional behavior when L scales as e w . 
Namely, if each bond is randomly occupied with proba- 
bility p, a percolation threshold at a value of p c strictly 
between and 1 arises, and associated critical exponents 
can be easily computed 0. 

To adapt this model to filtration, we posit that each 
bond has a radius r drawn from a specified distribu- 
tion, with volumetric flow rate proportional to r 4 Vp 
(Poiseuillc flow), where Vp is the pressure gradient in 
the bond. Dynamically neutral suspended particles move 
through the medium at a rate governed by this local flow. 
We assume perfect mixing at each node, in which a sus- 
pended particle has a flow induced probability propor- 
tional to rf to enter an unblocked bond of radius fj in 
the next downstream bundle. The particles, whose radii 
are also drawn from a distribution, are injected singly 
and tracked until each is trapped or escapes the system. 
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Upon capture, the particle is defined to block the bond 
completely so that there is no further fluid flow in this 
bond. After each blockage event, the new flow field is 
computed to determine the trajectory of the next sus- 
pended particle. 

This system exhibits three regimes of behavior. For 
pores typically smaller than particles (subcritical), the 
particles get trapped almost immediately and rapid clog- 
ging ensues. Conversely, for pores typically larger than 
particles (supercritical), a steady state is eventually 
reached for a finite length system, in which the smallest 
pores are blocked and the suspension flows freely through 
the remaining unblockable pores. These cases can be 
viewed as corresponding to poor filter performance. At 
the boundary between these regimes is the critical case, 
where the particle and pore radius distributions overlap 
substantially. Here, particle trapping is gradual, with 
considerable penetration of the medium before clogging 
is reached. This may be viewed as efficient filtration be- 
cause of the large number of particles filtered before clog- 
ging and the relatively long filter lifetime. Thus both 
from practical and theoretical perspectives, the critical 
case is the most interesting. 

For simplicity and concreteness, consider a uniform 
distribution of both particle and bond radii in the range 
[a, b]. More general continuous distributions can be 
straightforwardly treated, but little new qualitative in- 
sight emerges. Let us first determine the spatial distribu- 
tion of trapped particles during filtration. The gradient 
nature of the trapping process implies that the number 
of blocked bonds in downstream bubbles remains small, 
even at the percolation threshold (see Fig. 2). We there- 
fore employ an "unperturbed" approximation in which 
the initial bond radius distribution is continued to be 
used during the clogging process. Within this approxi- 
mation and assuming Poiscuillc flow, the probability that 
a particle of radius r gets trapped in a bubble is 
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for large w. (Exact calculation over all configurations of 
bond radii shows that the above large-w form is asymp- 
totically correct for w > 5.) Consequently, the probabil- 
ity that a particle gets trapped in the n th bubble is 

P„ = (1-P < )"- 1 P < . (2) 

Averaging over the distribution of particle radii gives, 
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Depending on the lower cutoff a, there are two differ- 
ent asymptotic behaviors for this trapping probability. If 
a = 0, the integral reduces to the beta function |hJ 
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where T(n) the gamma function, and for large n, (P n ) ~ 
0.1836 . . .n~ 6/5 . Conversely, if a ^ 0, then the asymp- 
totic behavior of the integral in Eq. (^|) , which arises from 
the contribution near the lower limit, can be written ap- 
proximately as 
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FIG. 2. (a) Trapping probability {Pn) versus n at the per- 
colation threshold for a bubble model of width w = 50 for co- 
incident and uniform distributions of particle and bond radii. 
Shown are data, based on 10 4 configurations, for (a, b) = (0, 1) 
(o) and (a,b) = (0.7, 1.0) (A), with the latter divided by 10 
for visualization. The straight lines have slopes —6/5 and 
— 2 respectively, (b) Trapping probability, based on inject- 
ing 4 x 10 6 particles into an unperturbed 500 x 1000 square 
lattice whose axes are oriented at 45° with respect to the av- 
erage flow, for coincident, uniform particle and bond radius 
distributions on (0, 1). The straight line has slope —1.27. 

From the denominator in the second line of Eq. (Q), the 
crossover between the n~ 6 / 5 and n~ 2 behaviors occurs 
when v(b 5 — a 5 ) < a 5 , or equivalently, n > n* = (b/a) 5 . 
While the exponents 6/5 and 2 are specific to the uniform 
radius distribution and the flow-induced bond entrance 
probability, the existence of the power law is generic and 
requires only the overlap of the bond and particle radius 
distributions. For example, for the Hertz distribution of 
particle and bond radii, p(r) = 2re~ r , the result corre- 
sponding to Eq. (||) is (P n ) oc n -4 / 3 . 

Qualitatively similar behavior for (P n ) occurs in lat- 
tice networks. In the spirit of our unperturbed approx- 
imation and to obtain relatively extensive data, we fo- 
cus on the spatial distribution of the initially injected 
particle. Later particles exhibit nearly the same spatial 
distribution of trapping location, but much more time 
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is needed for computing this distribution, since the net- 
work permeability must be recalculated after each trap- 
ping event. For n > 10, the best fit power law to the 
data is (P n ) ~ n^ 1 , with /j, » 1.27 (Fig. 2(b)). 

In the supercritical regime (bonds larger than parti- 
cles), (P n ) exhibits near-critical behavior, except that 
some particles can escape from the system. Conversely, 
in the subcritical regime (bonds smaller than particles), 
Eq. (|) gives (P„) oc exp[-n(a 5 - A 5 )/(B 5 - A 5 )], where 
(a, b) and (A, B) are, respectively, the ranges of the par- 
ticle and bond radius distributions. As a — A — > 0, the 
decay length (B 5 — A 5 )/ (a 5 — A 5 ) diverges and power law 
behavior of (P n ) is recovered. It is in this sense that coin- 
cident bond and particle radius distributions corresponds 
to a critical phenomenon. 

Let us now determine the number of particles that need 
to be injected to reach the clogging (percolation) thresh- 
old. For the bubble model, this means that all bonds in a 
single bubble are blocked. Since (P n ) monotonically de- 
creases in n, the probability that all w bonds are blocked 
in the n th bubble is non-zero only for small n. (Numeri- 
cally, for w = 100, for example, the probability of block- 
ing in bubble n is approximately 78.2%, 15.8%, 4.44%, 
and 1.21% for n = 1, 2, 3, and 4, respectively.) In the 
following, we therefore work within the approximation 
that it is only the first bubble that clogs, and that both 
the particle and bond radius distributions are uniform on 
[0,1]. 

We first compute the number of particles that need to 
be injected into a bond of radius < r < 1 before it 
becomes blocked p!lj ]. For N particles, the probability 
that all have their radii in the range [0, r] is simply r N . 
This can be re-interpreted as the probability that the 
maximum radius among N particles lies between and 
r. That is r N = J Q r P N (r')dr', with P N (r) the probabil- 
ity density that the maximum radius equals r. Conse- 
quently, -P/v(r) = Nr N ~ x , and the average radius of this 
largest particle is 
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Inverting this relation shows that of the order of (1 — r) _1 
particles need to be injected before a particle of suffi- 
ciently large radius enters to block a bond of radius r. 

Consider now a single bubble of w > 1 bonds. The 
number of particles needed to block bonds whose radii 
are in the range [r, r + dr] is Wjzz. Consequently, the 
total number of particles needed to block the bubble is 
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Here we again use Eq. (^) to determine that the largest 
and smallest bond radii in the bubble are r max jw 
and r m i n w 1/w, respectively. The integral is dominated 
by the behavior at the upper limit and gives 



Notice that a naive determination of the percolation 
threshold from N c (Pi) = w gives, using Eq. (0), N c = 
6w. The logarithmic factor in Eq. (§]) arises from the 
widest bonds for which many particles need to be in- 
jected before blocking occurs. 
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FIG. 3. Number of particles injected at percolation (nor- 
malized by w) versus lnu>. The number of configurations is 
10 5 for w < 10 2 , 



10 4 for 10 2 < w < 10 3 , and 10 3 for w > 10 3 . 



Simulations on the bubble model indicate that this log- 
arithmic w dependence for ^V c is independent of the pre- 
cise form of the entrance probability for a particular bond 
and similar details. Thus if the system size increases 
isotropically, the percolation threshold p c = 1 — N c / Lw 
approaches 1, with Lw is the total number of bonds 
in the system. To obtain a threshold value less than 
unity requires exponential anisotropy in which L ~ lnw. 
This result is qualitatively robust with respect to differ- 
ent particle and bond radius distributions. For example, 
for the Hertz distribution, following analogous computa- 
tions to those just outlined gives N c oc w(\nw) 2 . For 
the square lattice, on the other hand, simulations indi- 
cate that N c is linearly proportional to the system width. 
This corresponds to a percolation threshold which scales 
as 1 — 1/L. Thus cither L should be constant, or an al- 
ternative relation between the bond and particle radius 
distributions may be appropriate to define criticality for 
a finite-dimensional network. 

Finally, consider the behavior of the permeability dur- 
ing filtration, for which the bubble model again provides 
a useful starting description. Due to the series geometry 
of the bubble model, the permeability k can be written 
as 



(9) 



N c oc wlnw. 



(8) 



where k n denotes the permeability of the n th bubble. As 
the filter becomes constricted, the number of trapped 
particles in the n th bubble is proportion to n _/i (Eq. (Q)), 
with [i detail dependent. Numerically, we find that the 
permeability of the n th bubble also scales as a power of 
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the number of unblocked bonds, k n oc w — A/n v . Here v 
is also detail dependent and A is proportional to the total 
number of particles trapped, with A — > w corresponding 
to clogging (up to logarithmic factors). Thus the inverse 
permeability is 
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This integral is approximately constant and gives 
~ L/w, except close to clogging. To estimate the 
integral in this limit, note that for n close to one, the 
integrand is dominated by the divergence in the denom- 
inator, while for n > (A/iv) 1 ^ the second term in the 
denominator can be neglected. Splitting up the integral 
according to this prescription gives 
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with n* = (A/w) 1/l/ . The first integral may be estimated 
by defining v — wn v j A and then treating the resulting 
slowly varying factor of v 1 ^ in the numerator as constant 
compared to the divergent factor l/(v — 1). We thereby 
obtain 
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where the correction term in k" 1 is manifestly posi- 
tive near the clogging threshold (A — > w from above). 
This crude estimate shows that the permeability of an 
isotropic system (L oc w) is essentially unaffected by in- 
dividual bond blocking events until one bubble is nearly 
completely blocked, after which k discontinuously drops 
to zero. If, however, L ~ law, then k vanishes as a log- 
arithm in A — w, where A — w can be identified as the 
distance to the percolation threshold. 

In summary, a quasi-one-dimcnsional bubble model 
successfully describes various geometrical aspects of 
depth filtration. For coincident bond and particle ra- 
dius distributions, the number of particles trapped a dis- 
tance n downstream asymptotically varies as n _M , with 
/i dependent on details of these distributions. The per- 
colation threshold can be determined from extreme value 
considerations, and within the bubble model, the length 
must scale as the logarithm of the system width to have 
a percolation threshold strictly less than unity. For such 
a geometry, the permeability exhibits logarithmic depen- 
dence on (p — p c ) over a restricted range. These results 
may help explain previous simulations on relatively small 
systems [Q , where a threshold value close to unity and a 
permeability which varied rapidly in concentration near 
the threshold was observed. 

These unusual results are a consequence of the gradient 
nature of the particle trapping process which predomi- 
nantly affects the upstream portion of the network. This 
gradient aspect also has ramifications for efficient filter 



design. A system whose width is much greater than its 
length is needed to give a percolation threshold less than 
unity. This geometry has the desirable feature that a fi- 
nite fraction of the medium actually traps particles. On 
the other hand, the radius-average probability that a par- 
ticle escapes a system of length L vanishes as L 1_AI for 
(Pn) ~ n~^. Thus a small escape rate and a percolation 
threshold strictly less than unity, which ostensibly are 
the desired properties of a depth filter, cannot simulta- 
neously be satisfied in a spatially homogeneous medium. 
This suggests that an optimal filter should have a longitu- 
dinal varying local permeability which effectively cancels 
the gradient nature of particle trapping process. 
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